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Abstract 


This  report  summarizes  the  results  on  the  project  obtained  during  the  two  year  period.  A 
new  special  puipose  laboratory  equipment  called  the  Scattered  Light  Polariscope  has  been  designed 
and  constructed.  The  Polariscope  system  consists  of  a  device  for  loading  the  specimen,  optical 
alignment  devices  for  precise  positioning  of  the  specimen  with  respect  to  the  incident  light  and  an 
imaging  system  with  a  video  camera,  digitizing  frame  grabber  and  image  processing  software.  The 
systems  have  all  been  assembled  and  tested  for  proper  operation.  Certain  limitations  of  the  current 
photoelastic  equations  were  determined  and  therefore  a  complete  reexamination  of  the  photoelastic 
equations  was  performed.  Calibration  experiments  have  been  performed  to  determine  the  reliability 
of  the  technique  as  well  as  the  data  interpretation  procedures.  Results  have  been  obtained  on  the 
nature  of  the  stress  variation  near  the  tip  of  a  rounded  notch  and  a  crack.  Significant  differences 
have  been  demonstrated  between  the  two  and  three  dimensional  stress  fields. 
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INTRODUCTION 

The  stress  field  near  the  tip  of  a  through  crack  in  a  plate  of  finite  thickness  is  known  to  have  a  three 
dimensional  character.  However,  due  to  the  analytical  complexity  of  the  problem,  two  dimensional 
approximations  of  plane  stress  or  plane  strain  are  usually  introduced  in  interpreting  the  crack  tip 
stress  state.  Under  the  two  dimensional  approximation,  the  stress  field  near  the  crack  tip  is 
completely  determined  by  a  single  parameter  Kj,  called  the  stress  intensity  factor.  The  fracture 
toughness  of  a  material  is  then  that  value  of  Kj  at  which  crack  growth  initiates,  denoted  by  KIC. 
However,  crack  growth  does  not  occur  uniformly  through  the  plate  thickness  and  hence  the  two- 
dimensional  approach  represents  an  idealization  of  the  actual  situation.  Also,  the  critical  stress 
intensity  factor  at  initiation  seems  to  be  dependent  on  the  plate  thickness. 

An  attempt  was  made  by  Roaskis  and  Ravi-Chandar  (1984)  to  determine  experimentally, 
the  range  of  validity  of  the  two  dimensional  approximations.  Using  the  experimental  method  of 
caustics,  measurements  were  made  at  different  distances  from  the  crack  tip  and  compared  with  two 
dimensional  analyses.  It  was  found  that  a  state  of  plane  stress  was  established  only  at  distances  of 
about  one  half  of  the  plate  thickness  away  from  the  crack  tip.  Furthermore,  no  significant  region 
where  a  plane  strain  field  dominated  was  found.  Some  of  these  conclusions  have  also  been 
obtained  analytically  by  Yang  and  Freund  (1985).  The  main  upshot  of  these  results  is  that  the 
plane  asymptotic  singular  field  characterized  by  the  stress  intensity  factor  Kj  does  not  determine  the 
conditions  local  to  the  crack  tip  where  the  fracture  processes  occur.  It  is  still  true,  however,  that  Kj 
remains  a  useful  measure  of  the  energy  available  for  the  fracture  process  as  along  as  a  strain  energy 
density  function  is  used  to  characterize  the  material  behavior. 

Analyses  of  the  linear  elastic  problem  using  the  full  three  dimensional  theory  of  elasticity 
have  been  presented  by  Folias  (1975)  for  a  plate  of  finite  thickness  and  by  Bentham  (1977)  for  a 
quarter-infinite  crack  in  a  half-space.  However,  these  solutions  have  not  been  used  much  in 
fracture  mechanics,  probably  due  to  the  difficulties  in  interpretation.  Smith  and  co-workers 
(1975,1983,  and  1984),  have  used  frozen  stress  photoelasticity  to  investigate  the  three  dimensional 
problems.  While  much  of  their  early  work  concentrated  on  three  dimensional  crack  geometries 
such  as  the  part  through  and  comer  cracks,  recently,  they  have  addressed  inherently  three 
dimensional  problems  in  finite  thickness  plates.  Their  experimental  results  support  the  analysis  of 
Bentham  by  indicating  that  the  inverse  square  root  singularity  of  the  two  dimensional  analysis  is 
lost  near  the  faces  of  the  plate.  The  technique  of  frozen  stress  photoelasticity  involves  machining 
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of  the  specimens  into  thin  slices  after  stress  freezing  and  then  examining  the  slices  using  two 
dimensional  or  transmitted  light  photoelasticity.  In  the  present  investigation,  scattered  light 
photoelasticity  has  been  used.  This  is  a  non-invasive,  non -destructive  technique,  and  in  the 
present  investigation,  some  fundamental  improvements  have  been  made  in  both  the  interpretation 
of  the  data  and  the  implementation  of  the  experimental  scheme.  These  are  outlined  in  detail  in  the 
following  sections. 

II.  THE  THEORY  OF  PHOTOELASTICITY 

In  1939,  Weller  made  a  remarkable  discovery:  the  light  field  scattered  by  a  transparent  elastic 
material  subjected  to  a  stress  field  exhibited  fringes  that  could  be  related  to  the  stresses.  He 
published  this  information  as  a  half  page  letter  to  the  editor!  Drucker  and  Mindlin  (1940) 
developed  the  corresponding  theory  to  analyse  the  fringe  patterns  and  extract  the  stress  variation. 
Weller  (1941)  implemented  the  technique  of  scattered  light  photoelasticity  and  published  pictures  of 
fringes  that  were  of  astounding  quality;  very  few  pictures  of  similar  quality  appear  in  later 
literature.  Since  this  early  development,  a  number  of  authors,  notably,  Jessop  (1951),  Frocht  and 
Srinath  (1958),  Srinath  and  Frocht  (1962),  Cheng(1969),  and  Srinath(1983),  have  explored  the 
technique  further  both  in  theoretical  development  and  experimental  implementation.  However,  the 
technique  has  remained  esoteric,  not  only  because  of  specialized  equipment  and  interpretation 
needed  to  obtain  information,  but  also  because  the  sensitivity  (stress  fringe  value)  was  quite  low. 
In  the  following,  we  propose  a  novel  implementation  of  the  technique  of  scattered  light 
photoelasticity  that  enhances  the  sensitivity  considerably.  In  fact,  a  fringe  pattern  is  already 
established  in  the  specimen  and  applied  stresses  only  cause  changes  in  this  fringe  pattern,  thus 
making  the  technique  similar  to  mismatch  moire,  for  example.  We  briefly  re'  iew  some  aspects  of 
crystal  optics  in  Section  H.1  The  stress  optic  law  and  its  tensorial  character  are  discussed  in  detail 
in  Section  II.2.  The  principles  of  scattered  light  photoelasticity  and  its  inteipretation  are  outlined  in 
Section  n.3.  An  example  applied  to  the  proposed  model  material,  Homalite-100  is  given;  it  is 
demonstrated  that  this  material  is  an  optical  uniaxial  crystal.  Finally  in  Section  n.4,  the  new  model 
material  is  used  to  examine  the  classical  disk  problem  to  illustrate  the  application  and  then  the  notch 
and  crack  problems  are  considered. 
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II.  1  Some  Aspects  Of  Crystal  Optics 

Crystalline  media  exhibit  optical  anisotropy;  i.e.,  the  speed  of  propagation  of  light  depends  on  the 
direction  of  propagation  and  the  plane  of  polarization.  It  can  be  shown  that,  (see  Coker  and  Filon 
(1957),  Bom  and  Wolf  (1965))  for  any  given  direction  of  propagation,  two  polarizations  of  the 
electric  vector  characterizing  the  electromagnetic  wave  are  possible  and  in  general  the  speeds  of  the 
two  components  are  different.  This  fact  forms  the  basis  of  all  photoelastic  techniques.  We  shall 
briefly  review  the  effects  of  optical  anisotropy  resulting  from  electrical  anisotropy  in  crystalline 
materials.  Magnetic  anisotropy  and  its  effect  on  optical  anisotropy  are  small  in  most  materials  and 
therefore  is  not  considered  here.  The  magnetic  permeability  will  be  considered  to  be  unity. 

Electrical  anisotropy  may  be  represented  through  the  relationship  between  the  electric  field 
vector,  E,  and  electric  displacement  vector,  D,  as  follows1: 


Di =  eijEj  (1) 

where  e;j  are  the  components  of  the  dielectric  tensor.  The  inverse  relation  can  be  written  as 

Ej =  aijDj  (2) 


where  a;j  are  the  components  of  the  the  tensor  a,  which  is  the  inverse  of  e.  The  refractive  index 
can  be  denoted  by  a  tensor  n  given  by  Va.  The  "square-root"  of  a  symmetric  tensor  of  order  two 
is  uniquely  defined  when  in  the  common  principal  frame,  the  components  of  n  are  taken  as  the 
square  root  of  the  principal  values  of  a  (see  Stephenson  (1980)}.  In  general,  the  electric 
displacement  is  not  in  the  same  direction  as  the  electric  field.  In  such  anisotropic  media,  is  possible 
to  propagate  "two  monochromatic  plane  waves,  with  two  different  linear  polarizations  and  two 
different  velocities"2  in  any  given  direction.  Furthermore,  the  components  of  the  dielectric 
tensor  may  be  represented  graphically  as  an  ellipsoid 

EijXiXj  =  0  (3) 

This  is  referred  to  as  the  ray  ellipsoid  or  Fresnel  ellipsoid.  A  similar  representation  can  be 
provided  for  the  components  of  the  tensor  a;  the  corresponding  ellipsoid  is  called  the  index 
ellipsoid  or  the  ellipsoid  of  wave  normals.  In  light  of  the  new  definition  of  the  refractive  index  as 
a  tensor  square  root  of  a,  a  refractive  index  ellipsoid  could  also  be  constructed.  The  Fresnel 
and  the  index  ellipsoids  provide  a  graphical  means  of  determining  the  refractive  indices  for  rays 
propagating  in  arbitrary  directions.  Consider  a  ray  propagating  in  the  direction3  s  as  indicated  in 


'  Latin  subscripts  admit  the  range  1,2,3;  Greek  subscripts  admit  the  range  1,2.  Summation  convention  is  used. 
2  Born  and  Wolf  (1965),  p.671. 
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Figure  1 ,  where  a  Fresnel  ellipsoid  is  also  illustrated.  The  intersection  of  the  plane  normal  to  s 
with  the  ellipsoid  is  an  ellipse;  let  the  directions  of  the  major  and  minor  axes  of  this  ellipse,  which 
will  be  referred  to  as  the  secondary  principal  directions ,  be  denoted  by  x'a.  The  two  plane 
polarized  waves  that  can  be  propagated  lie  along  the  secondary  principal  directions.  The  speeds  of 
propagation  for  the  two  different  polarizations  are  given  by  the  semi-major  and  semi-minor  axes  of 
the  ellipse.  The  refractive  index  is  given  directly  by  the  secondary  principal  values  of  the  refractive 
index  tensor  n. 

Since  e  is  symmetric,  there  exists  a  principal  frame  in  which  the  components  of  e  can  be 
represented  in  diagonal  form.  Let  £j  be  the  principal  values  of  E.  If  all  the  principal  values  are 
equal,  then  =  eS-,  a-  =  aS^, where  £=l/a  and  £  and  a  are  constant.  Then  the  electric  field  and 
displacement  are  in  the  same  direction  and  the  medium  is  optically  isotropic;  the  corresponding 
Fresnel  ellipsoid  collapses  to  a  sphere.  If  two  principal  values  are  equal,  then  isotropy  is  lost;  the 
corresponding  Fresnel  and  index  ellipsoids  are  ellipsoids  of  revolution.  Thus,  the  direction 
corresponding  to  the  unique  principal  value  is  fixed  while  the  other  two  principal  directions  may  be 
chosen  arbitrarily  in  the  plane  perpendicular  to  the  unique  principal  direction.  This  unique  principal 
direction  is  called  the  optic  axis ;  such  a  medium  is  referred  to  as  a  uniaxial  crystal.  We  shall 
show  later  that  in  photoelastic  applications,  an  initially  isotropic  medium  becomes  a  uniaxial 
crystalline  medium.  If  all  three  principal  values  are  distinct,  then  there  are  no  symmetries 
associated  with  the  Fresnel  ellipsoid.  However,  it  can  be  shown  that  there  exist  two  distinct  axes 
normal  to  which  the  ellipsoid  is  sectioned  into  circles.  These  two  axes  are  the  optical  axes  for  this 
medium  and  materials  which  exhibit  this  behavior  are  referred  to  as  biaxial  crystals. 

Some  polymeric  materials  that  are  optically  isotropic,  become  anisotropic  upon  an  imposed 
stress  or  deformation.  If  the  deformation  is  reversible,  the  anisotropy  disappears  on  the  removal  of 
the  load;  on  the  other  hand,  if  large  irreversible  deformations  are  imposed,  the  anisotropy  becomes 
permanent.  In  other  polymers,  during  the  casting  or  molding  process,  a  permanent  molecular 
orientation  is  introduced  which  results  in  optical  anisotropy.  This  temporarily  or  permanently 
induced  anisotropy  can  be  used  to  determine  the  stress  and  deformation  in  such  materials  or  the 
molecular  orientation  and  the  residual  stress  field  and  this  idea  forms  the  basis  of  photoelasticity. 
In  the  following  section,  we  describe  the  stress  optic  law  that  is  based  on  Wertheim's  model. 


3  The  ray  direction  will  always  be  denoted  by  s. 


11.2.  Stress-Optic  Law 

In  this  section,  we  derive  the  forms  of  the  stress  optic  law  commonly  used  in  photoelasticity  in 
such  a  way  that  it  can  be  applied  to  transmitted  as  well  as  scattered  light  photoelasticity.  Viewing 
the  photoelastic  medium  as  an  optical  medium  that  may  be  initially  optically  isotropic  or  anisotropic 
and  which  becomes  weakly  optically  anisotropic  upon  the  application  of  load,  photoelasticity  is 
described  in  the  setting  of  crystal  optics  discussed  above.  This  is  essential  since  in  our 
implementation  of  scattered  light  photoelasticity,  we  use  an  uniaxial  optical  crystal  as  the  model 
material. 

Consider  a  light  ray  propagating  along  a  direction  s  as  indicated  in  Figure  1,  entering  the 
photoelastic  medium  subjected  to  a  stress  field  a(x).  Under  conditions  of  small  deformations,  it 
may  be  assumed  that  the  changes  in  the  dielectric  tensor  or  the  index  tensor  are  a  linear  function  of 
the  stress  state  or  deformation.  The  coefficients  in  this  relationship  are  referred  to  as  the  stress- 
optical  or  piezo-optical  coefficients.  This  linear  relationship,  suggested  by  Wertheim(1854), 
Maxwell(1852)  and  Neumann  (1841)  can  be  written  in  a  generalized  form  as  follows: 

aij(x)  -  ay  =  Xijki°ki(x)  (4) 

where  x  is  a  constant  tensor  of  fourth  order  called  the  stress-optic  tensor  with  symmetries 
somewhat  similar  to  the  elasticity  tensor  C  relating  the  stress  and  strain  tensors.  The  symmetry  of 
a  and  a  impose  certain  symmetries  on  x,  reducing  the  number  of  constants  in  the  most  general  case 
to  36.  Bhagavantham  (1942, 1952)  has  investigated  the  effect  of  various  crystal  symmetries  on  the 
tensor  X-  We  will  look  at  only  the  special  case  of  isotropic  photoelastic  symmetry  in  the 
photoelastic  effect;  i.e.,  x  is  isotropic.  The  medium,  however,  may  be  initially  optically 
anisotropic  or  isotropic,  and  the  applied  stress  a(x)  may  also  be  anisotropic,  resulting  in  optical 
anisotropy  in  the  stressed  medium.  Given  that  the  changes  in  the  tensor  a  under  the  applied  stress 
is  small,  Eq.(4)  may  be  replaced  with  a  corresponding  expression  for  the  refractive  index  tensor: 

nij(x)  -  n°ij  =Vijki°kl(x)  (5) 

In  fact,  since  experiments  only  measure  the  optical  path  difference  which  is  proportional  to  the 
difference  in  refractive  index,  the  above  equation  represents  a  more  convenient  form  of  the  stress 
optic  law. 

11.2.1  Isotropic  Photoelastic  Effect :  For  a  photoelastically  isotropic  medium,  i.e.,  the  stress  optic 
effect  is  independent  of  direction  of  propagauon  of  light,  only  two  components  of  X  are  unique  and 
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can  be  represented  as 

Xijld  —  P(8ikSj. 

where  a  and  P  are  the  only  stress  optic  coefficients  needed  to  define  the  photoelastic  effect. 
Substituting  the  above  into  Eq.(4),  yields: 

ay (x)  -  ag  =  ocolck(x)8ij  +  Poy(x)  (7) 

This  relationship  is  similar  in  form  to  the  linear  elastic  relationship  between  the  components  of  the 
strain  tensor  and  the  components  of  the  stress  tensor.  The  form  of  this  equation  indicates  that  a 
determines  the  hydrostatic  contribution  while  p  provides  the  contribution  of  the  (deviatoric)  shear 
stresses. 

Initially  isotropic  medium  with  isotropic  photoelastic  effect:  From  Eq.(7),  it  is  clear  that  the 
principal  axes  of  the  stress  tensor  and  the  principal  axes  of  the  refractive  index  tensor  in  the 
stressed  medium  coincide  if  and  only  if  the  medium  is  initially  isotropic ;  i.e.,  a0;j  =  a0  8^.  Then 
the  above  relationship  reduces  to  the  standard  form  of  the  stress  optic  law  when  expressed  in  this 
common  principal  frame. 

aj(x)  -  a0  =  (a  +  pja^x)  +  a[o2(x)  +  o3(x)] 

a2(x)  -  a0  =  (a  +  p)a2(x)  +  a[o3(x)  +  O^x)]  (8) 

a3(x)  -  a0  =  (a  +  p)a3(x)  +  ato^x)  +  o2(x)] 

where  (a  +  2P)  and  a  are  called  the  direct  and  transverse  stress  optic  coefficients  respectively. 
For  the  ray  propagating  in  an  arbitrary  direction  s,  the  incident  polarized  ray  is  decomposed  into 
two  components  and  propagated  with  different  velocities.  The  appropriate  refractive  indices  can  be 
obtained  through  the  graphical  construction  outlined  earlier.  Consider  a  right  handed  coordinate 
system  given  by  (s,  x'j  and  x’2)  as  indicated  in  Figure  1,  where  x'j  and  x'2  are  chosen  such  that 
they  are  along  the  semi-major  and  semi-minor  axes  of  the  ellipse  given  by  the  intersection  of  the 
plane  normal  to  s  with  the  Fresnel  ellipsoid.  Only  light  rays  whose  D  vectors  are  aligned  along  x’, 
and  x'2  can  propagate  along  the  s  direction.  The  corresponding  components  of  a  in  the  frame  (s, 
x’j  and  x'2)  are  a'n  and  a’22;  by  the  special  choice  of  the  coordinate  frame,  a'12  =  0.  The 
components  of  the  stress  tensor  o(x)  in  this  frame  are  denoted  by  a'n,  a’22,  a'ss,  o'^  and  <j's1; 
once  again,  a'I2  =  0  due  to  the  choice  of  coordinates.  The  components  o'n,  c'^,  and  o'ss  are 
referred  to  as  the  secondary  principal  stresses  and  are  the  only  components  that  influence  the 
changes  in  the  refractive  index4.  The  stress  optic  law  can  then  be  written  in  in  terms  of  the 
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secondary  principal  stresses: 

a'„(x')  -  a0  =  (a  +  (J)a'„(x')  +  a[o'22(x')  +  o'ss(x')l 
a'22(x')  -  a0  =  (a  +  p)a'22(x')  +  a[a'n(x')  +  a'ss(x')] 


Thus,  for  an  arbitrary  direction  of  propagation,  s,  two  component  rays  with  polarizations  along  the 
directions  x’j  and  x’2  are  propagated.  The  refractive  indices  for  these  two  component  rays  are 
given  by  n’t  =  1/Va'n  and  n'2  =  1/Va'22.  Thus,  the  two  rays  accumulate  a  phase  difference  as  the 
ray  propagates  through  the  medium.  The  relative  difference  is  given  by 


a'n(xO  -  a'22(x')  =  p  [o'„(x')  -  a'22(x')] 


(10) 


Note  that  the  relative  difference  is  independent  of  the  hydrostatic  component  of  stress.  Thus,  by 
measurements  of  the  relative  changes,  the  full  stress  field  cannot  be  determined.  However,  if  the 
absolute  changes  in  the  refractive  index  are  measured,  this  limitation  of  photoelasticity  may  be 
overcome.  In  terms  of  the  refractive  index,  Eq.(10)  can  be  written  as 

n'jCx')  -  n'2(x')  =  C  [a'nfx')  -  a'22(x')]  ^  ^ 

where  n’jfx’)  and  n’2(x')  are  the  refractive  indices  for  the  two  polarizations  that  can  propagate 
along  s.  In  terms  of  optical  path  difference,  the  "relative  retardation,"  A(x),  can  be  written  as: 

X  X 

A(x)  =  tn'jCx')  -  n'2(x')3ds  =  C[o'n(x')  -  &22(x)]  ds  (12) 


where  X  is  the  wavelength  of  light.  It  is  assumed  at  this  stage  that  the  secondary  principal 
directions  do  not  change  along  the  path  of  the  ray,  although  this  assumption  can  be  relaxed  with 
considerably  more  effort.  The  key  idea  in  photoelasticity  is  to  combine  the  two  components  and 
extract  the  phase  difference  for  determination  of  the  secondary  principal  stresses.  This  can  be  done 
in  two  ways:  first,  the  after  the  two  components  exit  the  photoelastic  medium  they  are  combined  by 
passing  through  a  polarizing  filter,  resulting  in  constructive  or  destructive  interference  depending 

4  If  the  stress  state  at  the  point  of  entry  is  described  by  the  Lame  stress  ellipsoid,  OjjXjXj  =  0,  then  the 

intersection  of  the  Lame  ellipsoid  with  the  plane  normal  to  s  describes  an  ellipse.  The  semi-major  and  semi-minor 
axes  of  this  ellipse  are  referred  to  as  the  secondary  principal  axes,  denoted  by  x'a  and  the  stress  components  along 

these  directions  are  referred  to  as  the  secondary  principal  stresses  ( effective  stresses  in  the  terminology  of  Weller), 
denoted  by  o'a.  These  definitions  are  similar  to  the  definition  of  the  secondary  principal  refractive  indices  and  the 

corresponding  axes.  If  the  standard  Wertheim's  law  of  temporary  double  refraction  (isotropic  photoelastic  effect)  is 
used,  then  the  ellipsoid  of  refractive  index  is  aligned  parallel  to  the  Lame  ellipsoid  of  stress.  Since  the  refractive 
indices  for  a  ray  propagating  in  the  direction  is  given  by  the  semi-major  and  semi-minor  axes  of  the  ellipse 
given  by  the  intersection  of  the  index  ellipsoid  with  the  plane  x^  =  0  (see  Bom  and  Wolf  (1965),  p.674),  the 
refractive  indices  for  rays  polarized  along  the  secondary  principal  stress  axes  are  a  linear  function  of  the  secondary 
principal  stresses  o'a. 
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on  the  phase  difference.  This  is  the  standard  application  of  transmitted  light  photoelasticity.  The 
resulting  fringe  patterns  can  only  provide  information  the  integral  of  the  secondary  principal 
stresses  along  the  ray  path. 

n  -  d'22  =  2^A(xl5X2)  =  -^N(Xl,x2)  (13) 


where  ct'n,  o'22  are  averaged  over  the  path  of  the  light  beam.  In  the  second  implementation,  the 
two  components  are  made  to  scatter  from  the  specimen  in  a  direction  normal  to  s.  This  scattering 
may  be  achieved  either  by  molecular  scattering  or  by  introducing  secondary  particles  in  the  model 
material  that  enhances  scattering.  Upon  scattering,  the  two  components  are  plane  polarized  in  the 
X’]-X'2  plane  and  combine  to  provide  interference  fringes.  Since  scattering  can  occur  from 
arbitrary  points  in  the  specimen,  the  thickness  variation  of  the  stress  components  can  be  extracted 


from  the  gradient  of  the  fringe  variation: 


cr'nfxO  -  o'22(x')  =  fa 


8N(x’) 

8s 


fa  8A(x’) 
2jch  8s 


(14) 


Thus,  through  the  scattered  light  implementation  of  photoelasticity,  it  is  possible  to  obtain  the 
spatial  variation  of  a’n(x')-  a'22(x')  without  thickness  averaging. 

Note  that  as  outlined  here,  the  only  difference  between  scattered  light  photoelasticity  and 
transmitted  light  photoelasticity  occurs  in  the  way  in  which  the  two  polarized  components  are 
recombined  to  provide  information  on  the  path  difference  A(x').  In  the  scattered  light  technique, 
the  specimen  itself  acts  as  the  'analyser'  used  in  transmitted  light  photoelasticity,  permitting 
evaluation  of  A(x')  even  in  the  interior  of  the  specimen. 


Initially  anisotropic  medium  with  isotropic  photoelastic  effect:  Now  let  us  consider  a  medium  for 
which  a0-  *  a0  5-.  Then  the  principal  axes  of  the  Lame  ellipsoid  and  the  Fresnel  ellipsoid  do  not 
coincide!  Eqs.  (4)  cannot  be  simplified  further.  However,  if  we  consider  the  light  propagation 
direction  s  to  lie  completely  in  a  plane  of  symmetry  of  the  specimen,  then  one  principal  direction 
must  lie  along  s,  implying  that  and  c’sl  are  zero.  c'H,  a22,  o'ss  and  0’12  are  nonzero.  The 
secondary  principal  stresses  cfj,  a'2,  can  be  obtained  from  Mohr's  transformation  as 


and  the  stress  optic  law  becomes: 


9 


(16) 


a'n(x')  -  a?!  =  (a  +  f^U')  +  a[a'2(x')  +  a's(x')l 
a'22(x')  -  a22  =  (<x  +  p)o'2(x')  +  a [a'^x')  +  a's(x')l 
Once  again,  if  only  the  relative  retardation  is  considered,  we  get 

a'n(x')  -  a'22(x')  =  {a'?,  -  a^2}  +  (3{o'1(x/)  -  a'2(x')>  (17) 

In  terms  of  the  refractive  indices,  the  relative  retardation  may  be  written  as 

n',(x')  -  n'2(x')  =  (n^-  n^}  +  C{o/1(x')  -  a'2(x')}  (18) 

Thus  the  stress-optic  law  for  transmitted  and  scattered  light  photoelasticity  may  be  written  as: 

X  X 

A(x)  =  y-J  [n'jCxO  -  n'2(x')]ds  =  —{rT?  -  n^x  +  C[a'„(x')  -  o'22(x')]  ds 

(19) 

For  planes  other  than  those  of  symmetry,  Eqs.(4)  must  be  used  and  detailed  measurements  of  the 
refractive  index  have  to  be  made  in  order  to  identify  the  principal  axes  and  thereby  identify  the 
stress  components.  We  do  not  investigate  this  aspect  here,  but  will  address  it  in  a  future 
contribution.  We  now  turn  to  scattered  light  photoelasticity  and  its  implementation  in  an  initially 
anisotropic  medium. 

III.  SCATTERED  LIGHT  PHOTOELASTICITY:  A  NOVEL  IMPLEMENTATION 
The  details  of  the  implementation  of  this  technique  are  available  in  text  books  such  as  those  by 
Dally  and  Riley  (1978)  and  Kukse  and  Robertson  (1974).  A  more  recent  monograph  by  Srinath 
(1983)  is  devoted  entirely  to  scattered  light  photoelasticity.  Therefore,  we  do  not  provide  an 
analysis  of  the  light  vector  as  it  travels  through  the  specimen;  only  the  results  are  given.  Consider 
a  sheet  of  light  polarized  along  the  direction  1-1  travelling  in  the  direction  s  as  shown  in  Figure  2. 
Let  the  secondary  principal  direction  x'2  be  oriented  at  an  angle  0  with  respect  to  the  plane  of 
polarization  1-1.  If  the  light  scattered  from  the  point  Q  is  viewed  along  a  direction  at  an  angle  <p  to 
the  secondary  principal  direction  x'2,  then  the  intensity  of  the  scattered  light  field  may  be  obtained 
as: 

I(x')  =  K  {cos20sin2(p  +  sin^cos2^  -  (1/2)  sin20  sin29  cosA(x'))  (20) 

where  K  depends  on  the  incident  light  amplitude  as  well  as  the  scattering  properties  of  the  medium. 
It  is  easily  shown  that  when  0  =  tp  =  ti/4,  the  resulting  intensity  is  given  by 

I(x')  =  (1/2)K  (1  -  cosA(x'))  (21) 
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Thus,  if  the  secondary  principal  directions  are  known  (fixed),  then  the  angles  <p  and  0  may  be 
adjusted  such  that  the  scattered  light  field  intensity  depends  only  on  the  secondary  principal 
stresses.  The  fringes,  formed  as  the  loci  of  the  corresponding  light  intensity  points,  will  be  similar 
to  the  isochromatics  in  transmitted  light  photoelasticity,  although  the  interpretation  should  the 
through  Eq.  (19).  On  planes  of  symmetry  and  on  stress  free  bounding  surfaces,  the  secondary 
principal  directions  can  easily  be  established  and  therefore,  we  illustrate  in  the  present  paper  a 
applications  of  scattered  light  photoelasticity  under  these  conditions. 

III.  1  Scattered  Light  Polariscope 

The  conventional  transmitted  light  polariscope  is  not  suitable  for  the  scattered  light  investigations. 
Weller  (1939)  built  a  rather  simple,  inexpensive  version  of  the  scattered  light  polariscope, 
concentrating  mainly  on  the  generation  of  the  incident  light  sheet,  with  the  appropriate  plane  of 
polarization.  Most  polariscopes  built  later  did  not  use  light  sheets,  but  merely  light  pencils,  thereby 
limiting  the  available  information  from  any  one  experiment  Weller  also  pointed  out  to  the  need  for 
immersing  the  complete  specimen,  along  with  any  loading  frame,  in  a  tank  filled  with  a  fluid  of 
matching  refractive  index,  to  eliminate  refractions  at  entry  into  and  exit  from  the  specimen.  The 
polariscope  built  for  the  present  investigation  is  shown  in  Figure  3.  The  main  idea  used  in  the 
design  of  this  polariscope  is  that  the  plane  of  the  light  sheet  remains  fixed  in  space;  the  plane  of 
polarization,  the  plane  of  observation  and  the  orientation  of  the  specimen  which  is  intersected  by 
the  light  sheet  are  variable  and  controlled  by  precision  rotating  and  translating  stages.  A  6mW  He- 
Ne  Laser  is  used  as  the  coherent  light  source  which  produces  a  pencil  of  light  rays  with  a  known 
plane  of  polarization.  The  laser  beam,  initially  travelling  horizontally,  is  deflected  upwards  along 
the  direction  s,  by  a  mirror.  The  light  beam  is  then  expanded  and  collimated  into  a  light  sheet  by  a 
pair  of  cylindrical  lenses.  The  laser  beam  can  be  rotated  about  the  s  axis,  thereby  permitting 
variation  of  the  plane  of  polarization  with  respect  to  the  plane  of  the  light  sheet  The  specimens  are 
loaded  in  a  specially  designed  loading  frame,  containing  a  hermetically  sealed  load  cell,  that  can  be 
immersed  in  a  tank  containing  a  fluid  with  a  refractive  index  matching  that  of  the  specimen. 

1 1 1.2  Scattering  in  Optical  Crystals 

Consider  a  light  ray  propagating  along  the  direction  s  incident  on  an  optically  anisotropic  medium, 
linearly  polarized  at  an  angle  <p  with  respect  to  the  major  principal  direction.  Upon  entering  the 
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medium,  the  polarized  light  is  decomposed  into  two  components  each  polarized  along  the  principal 
directions.  Since  the  speed  of  propagation  of  the  two  components  are  different,  they  accumulate  a 
path  difference  A,  given  by  Eq.  (19)  with  a  =  0.  Upon  scattering,  the  intensity  is  given  by 


I(x')  =  —K\ 


1-cos  -  nfy{s  -  {)}  j 


(22) 


where  {s-^(x’j)}  =  0  represents  the  equation  for  the  entrance  boundary  of  the  specimen  in  the 
plane  x'2  =  0..  Thus  the  fringes  should  be  parallel  to  the  boundary  of  the  specimen.  Scattering 
properties  of  most  optical  crystals,  such  as  Nicols  prism,  is  very  poor  and  thus,  the  prediction  of 
Eq.  (22)  is  not  easily  observed.  However,  Homalite-1005,  a  thermosetting  polyester,  possesses 
good  molecular  scattering  properties  and  is  also  optically  anisotropic.  The  latter  is  due  to  molecular 
orientation  introduced  in  the  casting  process.  The  prediction  of  Eq.(22)  is  shown  to  be  correct 
experimentally,  as  illustrated  in  Figure  4,  where  the  scattered  light  fields  due  to  light  incidence  on 
three  different  planes  of  a  circular  disk  are  shown. 

If  the  scattered  light  intensity  predicted  by  Eq.(22)  is  measured,  then  such  measurements 
can  be  used  to  determine  the  components  of  the  refractive  index  tensor.  This  is  illustrated  below. 


III. 3  Characterization  of  Optical  Properties  of  Homalite-100 

We  now  turn  to  an  investigation  of  the  optical  anisotropy  of  Homalite-100,  in  order  to  (i)  determine 
the  principal  refractive  indices  and  (ii)  to  determine  the  stress-optic  constants  through  some 
calibration  experiments.  Since  the  specimen  material  is  manufactured  in  the  form  of  cast  sheets, 
we  shall  introduce  a  consistent  coordinate  notation  to  keep  track  of  the  sheet,  the  principal 
refractive  axes  and  the  principal  stress  axes.  The  direction  normal  to  the  sheet  of  the  specimen  will 
be  denoted  by  X3.  The  directions  in  the  plane  of  the  specimen  will  be  denoted  by  Xa  and  the 
specific  orientation  will  be  dictated  by  the  specimen  boundaries  in  the  plane  X3  =  0.  The  principal 
refractive  axes  will  be  denoted  by  Nj,  with  the  corresponding  refractive  indices  indicated  by  n;. 
The  principal  stress  axes  will  be  denoted  by  Zj  and  the  corresponding  stress  components  denoted 
by  o ;.  The  coordinate  system  specific  to  the  light  propagation  direction  will  still  be  indicated  by  (s, 
x'j  and  x’2)  where  s  is  the  light  propagation  direction  and  x\  and  x'2  are  the  secondary  principal 
directions. 

III. 3.1.  Optical  Anisotropy  of  Homalite-100:  Figure  4  shows  three  photographs  of  the  scattered 
light  field  under  the  following  conditions:  (i)  s  direction  is  aligned  parallel  to  Xj  with  the  light 

5  Trademark  of  SGL  Homalite,  Delaware 
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sheet  in  the  X,-X2  (4a)  and  X3-X,  (4b)  planes;  (ii)  s  direction  is  aligned  parallel  to  X3  with  the 
light  sheet  in  the  X3-Xj  (4c)  plane.  First,  note  that  the  existence  of  variations  in  the  intensity  of  the 
scattered  light  field  in  the  different  pictures  is  indicative  of  the  optical  anisotropy  of  Homalite-100. 
Secondly,  when  s  is  in  the  direction  of  X3,  uniform  scattered  light  intensity  is  observed.  This 
implies  that  the  rays  with  polarizations  along  the  secondary  principal  directions  x'j  and  x'2  do  not 
experience  any  path  difference;  i.e.,  the  secondary  principal  refractive  indices  n'j  and  n'2  are 
equal.  This  is  possible  in  a  general  anisotropic  medium  (biaxial  crystal)  only  if  s  were  in  the 
direction  of  the  optic  axis.  Finally,  by  continuously  varying  the  direction  s  from  the  X3  direction 
to  -X  direction  in  the  X2-X3  plane,  no  other  optic  axis  was  found  implying  that  the  specimen  was 
a  uniaxial  crystal,  with  the  principal  refractive  index  axes  N;  aligned  along  the  specimen  axes  X;. 
i.e.,  the  index  ellipsoid  is  an  ellipsoid  of  revolution,  with  the  axis  of  revolution  aligned  along  X3  or 
N3. 

III. 3. 2.  Stress  Optic  Coefficients  of  Homalite-100:  The  stress-optic  coefficients  were  determined 
through  a  number  of  calibration  experiments  using  simple  specimen  loading  conditions  where  the 
stress  field  is  completely  known  from  elementary  elastic  analysis.  Planar  tensile  specimen  and  a 
disk  with  diametral  compression  were  used  in  the  calibration  experiments.  In  both  these  cases, 
plane  stress  conditions  may  be  assumed  since  the  plate  thickness  was  quite  small.  Thus,  the 
principal  stress  direction  Z3  coincides  with  the  direction  X3  of  the  specimen  and  also  with  the 
principal  refractive  index  direction  N3.  Furthermore,  the  uniaxial  crystalline  nature  of  the  specimen 
implies  that  when  a  plane  stress  field  is  applied,  the  principal  stress  axes  Zj  and  Z2  should 
coincide  with  the  principal  refractive  index  axes  Nj  and  N2  (planar  isotropy).  Thus  in  this  instance 
(and  in  any  case  when  there  is  a  plane  of  symmetry  that  is  aligned  with  a  principal  refractive  index 
axis)  the  Lame  ellipsoid  and  the  Fresnel  ellipsoid  are  parallel.  Two  stress-optic  constants 
indicating  the  changes  in  n'j  and  n'2  should  completely  determine  the  stress  optic  effect  for  s 
coinciding  with  X3.  Furthermore,  if  one  is  interested  only  in  the  relative  retardation  as  in  Eq.(19), 
then  only  one  stress  optic  coefficient  has  to  be  determined. 

Consider  a  bar  loaded  in  tension  as  indicated  in  Figure  5.  The  orientation  of  the  specimen 
axes,  loading  direction  and  the  principal  refractive  axes  are  also  indicated  in  this  figure.  The 
fringes  in  the  unloaded  state  are  parallel  to  the  entry  surface  and  are  uniformly  spaced.  From  the 
fringe  spacing,  once  again  n'j  -  n’3  can  be  determined.  Upon  applying  a  uniform  tensile  stress 
field,  the  fringe  spacing  changes,  but  still  remains  straight.  From  the  new  fringe  spacing,  using 
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Eq.(19),  the  applied  stress  o'22  can  he  obtained,  if  fc  is  known;  otherwise,  this  experiment  is  used 
to  calibrate  fc. 

IV.  APPLICATIONS  OF  SCATTERED  LIGHT  PHOTOELASTICITY 
We  now  turn  to  some  examples  to  illustrate  the  capability  of  the  technique  of  scattered  light 
photoelasticty.  Two  examples  will  be  illustrated;  the  first  one  is  the  classical  problem  of  a  circular 
disk  under  diametral  compressive  loads.  It  is  shown  that  by  observing  the  scattered  light  field 
from  three  orientations,  the  stress  components  oQp  can  be  determined.  Furthermore,  if  it  is 
assumed  that  plane  conditions  hold,  which  is  reasonable  for  this  geometry,  the  stress  components 
oap  can  be  obtained  directly  at  all  field  points  with  the  need  for  elaborate  separation  techniques 
used  in  transmitted  light  photoelasticity!  The  second  example  concerns  a  plate  with  a  notch  (similar 
to  a  compact  tension  specimen,  except  that  the  crack  is  replaced  with  a  rounded  notch).  Finally, 
the  stress  field  near  the  tip  of  a  crack,  and  its  dependence  on  the  specimen  thickness  is  investigated. 

IV.  1.  Determination  of  Stress  Components  from  Scattered  Light  Field 
We  now  turn  to  the  determination  of  the  individual  stress  components  from  observations  of  the 
scattered  light  field.  Consider  a  plane  of  symmetry  or  the  bounding  surface  of  a  plate,  with  some 
stress  distribution.  Clearly  o3ot  =  0;  on  the  bounding  surface,  o33  is  also  zero.  If  the  scattered 
light  from  a  light  sheet  lying  on  this  plane  is  observed,  the  intensity  is  given  by  Eq.(21)  where 
D(x)  is  given  by  Eq.  (19).  The  secondary  principal  stresses  will  depend  on  the  direction  of  light 
propagation  s  in  the  Xl  -  X2  plane.  Consider  three  angles  of  light  incidence  along  0  =  0  and  0  = 
±  a,  as  indicated  in  Figure  6.  The  corresponding  secondary  principal  stresses  are  *  033,  G'n  - 
o33  and  o"n  -  o33.  But,  a'n  and  o"n  are  related  to  through  the  Mohr  transformation 
equations: 

o'n  =  On  cos2a  +  o22  sin2a  +  2c12sina  cosa 

.  .  (23) 

o"n  =  Oji  cos2a  +  o22  sinza -  20i2sina  cosa 

Thus,  using  Eq.(14),  from  the  measurements,  the  fringe  gradients  corresponding  to  the  secondary 
principal  stress  differences  c22  -  o33,  o’n  -  o33  and  a"jj  -  o33  can  be  obtained.  Let  us  denote  the 
corresponding  fringe  gradients  as  follows: 
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.  ,  ,  ,  3N(x)  3N0 

Oj2(x)-o„(x)  =  -5j---5^-sG(x) 

dN'(x)  dN0 
O' 

<J"„(x)  -  033(x)  =  -  J£r-  =  G"(x) 

Then  the  stress  differences  are  obtained  as: 

022(x)  -  o33(x)  =  G(x) 

G'(x)  +  G"(x)  -2G(x) 


an(x)-a22(x)  = 


2  cos2a 


Gi2(x)  = 


G'(x)  -  G"(x) 
2  sin  2a 


(24) 


(25) 


once  the  variation  of  the  fringe  gradients  at  every  point  is  known.  If  the  plane  is  on  a  bounding 
traction  free  surface,  o33  is  zero  and  hence  the  above  equations  can  be  used  to  determine  the 
individual  stress  components  directly  without  the  need  for  elaborate  separations  techniques  used  in 
transmitted  light  photoelasticity. 


IV.  1.  Disk  under  Diametral  Compression 

As  an  example  of  the  application  of  the  technique,  consider  the  disk  problem  illustrated  in  Figure  7. 
The  unloaded  and  loaded  scattered  light  fields  for  once  angle  of  incidence  along  the  line  of  the  load 
application  is  shown.  If  the  fringe  gradients  before  and  after  loading,  along  the  light  path,  are 
determined,  then,  the  variation  of  the  stress  component  can  be  obtained.  Figure  8a  shows  the 
fringe  spacing  data  before  and  after  loading,  and  Figure  8b  shows  the  calculated  values  of  an;  for 
comparison,  the  classical  elastic  results  (Timoshenko  and  Goodier  (1970))  are  also  shown  in  this 
figure.  As  can  be  seen,  there  is  good  agreement  between  the  theory  and  experiment  It  is  possible 
to  determine  the  stress  components  at  every  point  completely  using  this  method,  but  the  above 
suffices  to  illustrate  the  point. 


IV. 2.  Notch  Specimen 

The  second  example  of  application  concerns  a  notched  specimen  shown  in  Figure  9.  Two  different 
specimen  thicknesses  were  considered  (9.525mm  and  12.7mm).  The  scattered  light  field  from  the 
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unloaded  and  loaded  specimen  are  shown  in  Figure  10,  for  light  incidence  along  the  -Xj  direction. 
From  these  and  two  other  sets  of  observations  of  the  scattered  light  fields  corresponding  to  0  =  ± 
a,  Eq.(25)  could  be  used  to  determine  on  -  a22  and  o12.  The  experiment  can  be  repeated  for  light 
incidence  along  both  the  specimen  midplane  and  the  free  surfaces.  Figure  1 1  shows  the  variation 
of  these  stress  components  along  the  xt  axis.  It  is  seen  that  the  mid  plane  stresses  are  considerably 
higher  than  the  face  plane  stresses.  In  Figure  12,  isochromatic  contours,  corresponding  to  equal 
values  of  Oj  j  -  c22  are  shown.  Figure  12a  shows  the  contour  corresponding  to  the  same  stress 
difference,  while  Figure  12b  shows  the  two  different  contours  correpsonding  to  two  different 
values  of  the  stress  difference  that  pass  through  the  same  point  ahead  of  the  notch  tip.  Further 
analysis  of  the  stress  fields  is  possible  using  the  technique  illustrated  here  to  determine  the  details 
of  the  three  dimensional  variation,  but  this  is  not  attempted  here. 

IV. 3.  Compact  Tension  Specimen 

We  consider  next,  a  Compact  Tension  (CT)  specimen,  in  order  to  investigate  the  three 
dimensionality  of  the  crack  tip  stress  field.  From  the  scattered  light  analysis,  the  stress  component 
g22  is  determined  directly  at  all  points  on  the  free  surfaces,  and  -  o33  is  determined  at  all  points 
on  the  specimen  mid-plane.  With  the  observation  technique  discusses  earlier,  it  is  not  possible  to 
separate  a22  and  a33  in  the  midplane  without  investigation  of  the  stress  field  in  other  planes.  The 
variations  of  g22  -  o33  on  the  midplane  and  G22  on  the  face  plane  with  distance  from  the  crack  tip 
are  shown  in  Figure  13.  In  Figure  14,  the  isochromatics  corresponding  to  (an  -  o22)  from  the 
midplane  and  face  plane  are  shown  for  the  two  specimen  thicknesses  of  6.35mm  and  12.7mm,  for 
the  same  value  of  (gh  -  o22).  Since  the  distance  at  which  this  ischomatic  exists  in  the  thicker 
specimen  is  much  closer  to  the  crack  tip  than  in  the  thinner  specimen,  we  see  that  Figure  13a, 
corresponding  to  the  6.35  mm  thick  specimen  shows  close  agreement  with  the  2D  theory,  while  te 
conrours  in  Figure  13b,  corresponding  to  the  12.7mm  thick  specimen,  exhibits  large  differences. 
From  an  examination  of  the  results  in  Figures  13  and  14,  it  is  immediately  apparent  that  the 
midplane  stresses  are  higher  than  the  face  plane.  Furthermore,  at  a  distance  of  r/h  =  0.5,  the  mid 
plane  and  face  plane  results  appear  close  to  the  2D  theory.  Thus,  the  range  of  validity  of  the  2D 
field  is  definitely  at  distances  greater  than  r/h  =0.5.  This  reinforces  earlier,  indirect  results  of 
Rosakis  and  Ravi-Chandar  (1986).  Moreover,  recent  numerical  computations  of  Parsons,  Hall 
and  Rosakis  (1986),  indicate  that  the  a22  near  the  face  plane  has  exactly  the  same  variation  shown 
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in  Figure  13b.  Their  results  are  reproduced  in  Figure  15.  Figure  15a  corresponds  to  the  specimen 
midplane  and  Figure  15b  shows  the  radial  variation  of  o22  in  the  traction  free  surface.  Their 
results  for  the  midplane  are  not  directly  comparable  since  they  plot  only  the  radial  variation  of  o22 
and  not  (a, ,  -  a22)  as  we  have  done  in  Figure  13a.  From  the  results  of  Parsons  et  al  (1988),  o33 
is  a  positive  quantity  which  decays  away  from  the  crack  tip  and  goes  to  zero  at  about  r/h  =  0.5;  also 
near  the  crack  tip,  as  once  proceeds  from  the  traction  free  surfaces  to  the  midplane,  a33  increases. 
Therefore,  we  expect  that  c22  in  the  midplane  will  be  larger  than  the  2D  prediction.  Further  work 
is  under  progress  to  determine  the  full  variation  of  the  stress  components  in  other  planes,  and  will 
be  reported  elsewhere. 
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Figure  1.  Fresnel  ellopsoid  representation  of  the  dielectric  tensor  in  an 
anisotropic  material. 


3.  Scattered  Light  Polariscope. 
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Light  direction 


Figure  5.  Scattering  arrangement  for  a  simple  tension  specimen. 
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Figure  8.  (a)  Fringe  number  variation  along  the  light  path  in  the  disk  for  the 
loaded  and  unloaded  specimen,  (b)  Secondary  principal  stress 
difference  calculated  from  the  fringe  patterns  compared  with  the 
two  dimensional  elasticity  solution. 


Figure  10.  Scattered  light  field  in  the  unloaded  and  loaded  notch  speci 
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Figure  1 1.  Variation  of  the  secondary  principal  stresses  in  the  specimen  midplane 
and  traction  free  surface,  for  two  specimen  thicknesses. 
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Figure  12. 
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(a)  Isochromatics  on  the  midplane  and  face  plane  for  the  same  value 
of  stress  difference,  (b)  Isochromatics  that  pass  through  the  same 
point  ahead  of  the  notch  tip  (the  corresponding  stress  differences  are 
marked  in  the  figure). 
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Figure  13.  Variation  of  the  normal  stress  ahead  of  the  crack  tip  (a)  in  the  specimen 
midplane  and  (b)  in  the  traction  free  surface. 
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Figure  14.  Isochromatics  corresponding  to  the  same  value  of  stress  difference 
(a)  on  the  face  plane  and  (b)  on  the  specimen  midplane 
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Figure  15.  Variation  of  the  normal  stress  ahead  of  the  crack  tip  in 

(a)  the  midplane  and  (b)  the  traction  free  surface  obtained 
by  Parsons  et  al. 


